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We analyze the Zeno phenomenon in quantum field theory. The decay of an un- 
stable system can be modified by changing the time interval between successive 
measurements (or by varying the coupling to an external system that plays the role 
of measuring apparatus) . We speak of quantum Zeno effect if the decay is slowed 
and of inverse quantum Zeno (or Heraclitus) effect if it is accelerated. The analysis 
of the transition between these two regimes requires close scrutiny of the features 
of the interaction Hamiltonian. We look in detail at quantum field theoretical 
models of the Lee type. 



1 Introduction 

The seminal formulation of the quantum Zeno effect by Misra and SudarshanS 
deals with unstable systems, i.e. systems that decay following an approxi- 
mately exponential law.B Such a formulation was implicit also in previous 
work,0 where the features of the "nondecay" amplitude and probability were 
investigated. The attention was diverted to oscillating systems, character^ 
ized by a finite Poincare time, when Cook published a remarkable paperp 
proposing to test the quantum Zeno effect (QZE) on a two-level system un- 
dergoing Rabi oscillations. Although oscillating systems are somewhat less 
interesting in this context, they are also much simpler to .analyze and indeed 
motivated interesting experimentsB and lively discussions ,u giving rise to new 
ideasE However, interesting new phenomena-occur when one considers unsta- 
ble systems, whose Poincare time is infinite JjoEj in particular, other regimes 
become possible, in which measurement accelerates liajd^mamical evolution, 
giving rise to an inverse quantum Zeno effect (IZE).ElrE3li3'Eil 

The study of Zeno effects for bona fide unstable systems requires the use of 
quantum field .theoretical techniques and in particular the Weisskopf-Wigner 
approximatioriij and the Fermi "golden" rule:Ej for an unstable system the 
form factors of the interaction play a fundamental role and determine the 
occurrence of a Zeno or an inverse Zeno regime, depending of the physical 
parameters describing the system. The occurrence of new regimes is relevant 



from an experimental perspective, in view of the beautiful experiments re- 
cently performed by Raizen's group on the short-time non-exponential decay 
(leakage through a confining potential)t3 and on the Zeno effects for such 
(nonoscillating) systemsJla 

We analyze here the transition from Zeno to inverse Zeno in a quantum 
field theoretical context, by looking in particular at the Lee model. This 
is a good prototype forxither field-theoretical examples and is more general 
than one might think.EjlHJ The usual approach to QZE and IZE makes use of 
"pulsed" observations of the quantum state. However, one obtains essentially 
the same effects by performing a "continuous" observation of the quantum 
state, e.g. by means of an intense field, that .plays the role of external, "mea- 
suring" apparatus. IJus is not a new idea,aa but has been put on a firmer 
basis only recently.Eil The 'icontinuous" formulation of the QZE has been 
discussed in detail elsewheretBEj and will be briefly reconsidered here, by fo- 
cusing in particular on quantum field theory and its interesting peculiarities, 
leading to new effects. 

2 Notation and preliminary notions: pulsed measurements 

Let H be the total Hamiltonian of a quantum system and |o) its initial state 
at t = 0. The survival probability in state \a) is 

P(t) = \A(t)\ 2 = \{a\e- iHt \a)\ 2 (1) 

and a short-time expansion yields a quadratic behavior 

P(t) ~ 1 - t 2 /rl r z 2 = (a\H 2 \a) (a\H\a) 2 , (2) 

where rz is the Zeno time. Observe that if the Hamiltonian is divided into a 
free and an interaction parts 

H = H +H int , with H \a) =u a \a), (a\H int \a) = 0, (3) 

the Zeno time reads 

r z 2 = (a\H 2 nt \a) (4) 

and depends only on the (off-diagonal) interaction Hamiltonian. _ 
Let us start from "pulsed" measurements, as in the seminal approach!!! 
The notion of "continuous" measurement will be discussed later. Perform N 
(instantaneous) measurements at time intervals r = t/N, in order to check 
whether the system is still in state \a). The survival probability after the 
measurements reads 

P(")(t) - P(t) n = P(t/N) N ~ exp H 2 /r|iV) N -^° 1. (5) 
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Figure 1. Evolution with frequent "pulsed" measurements: quantum Zeno effect. The 
dashed (full) line is the survival probability without (with) measurements. The gray line is 
the interpolating exponential (ph. 



If N = oo the evolution is completely hindered. For very large (but finite) N 
the evolution is slowed down: indeed, the survival probability ^iter N pulsed 
measurements (t = Nt) is interpolated by an exponential lawE3 

pW(i) = P(r) N = exp(iV log P(r)) = exp(- 7eff (r)t), (6) 

with an effective decay rate 

7eff(r) = -ilogP(r) = --log|^(r)| = -~Re[log A(r)} > . (7) 

T T T 

For r — > one gets P(t) ~ exp(— t 2 /t|), whence 

7cff(r) ~ r/r| (r 0) (8) 

is a linear function of r. Increasingly frequent measurements tend to hinder 
the evolution. ThxL physical meaning of the mathematical expression "r — > 0" 
is a subtle issue,E3c3 that will be touched upon also in the present article. 
The Zeno evolution is represented in Figure |l|. 



3 From quantum Zeno to inverse quantum Zeno ("Heraclitus") 

Let us concentrate our attention on truly unstable, systems, with a "natural" 
decay rate 7, given by the Fermi "golden" rule.E3 We ask: is it possible to 
find a finite time r* such that 

7eff(r*)= 7 ? (9) 

If such a time exists, then by performing measurements at time intervals t* 
the system decays according to its "natural" lifetime, as if no measurements 
were performed. The quantity r* is related to Schulman's "jump" time.Ej'B 
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Figure 2. Z < 1. (a) Determination of r* . (b) Quantum Zeno vs inverse quantum Zeno 
( "Heraclitus" ) effect. 



By Eqs. (|) and (|) one gets 

P(r*) = e-T T *, (10) 

i.e., t* is the intersection between the curves P(t) and e _7t . In the situation 
outlined in Figure ||(a) such a time t* exists: the full line is the survival 
probability and the dashed line the exponential e -7 * [the dotted line is the 
asymptotic exponential Ze~ 7 ', that will be defined in Eq. (plj)]. The phys- 
ical meaning of r* can be understood by looking at Figure Wb), where the 
dashed line represents a typical behavior of the survival probability P(t) when 
no measurement is performed: the short-time Zeno region is followed by an 
approximately exponential decay with a natural decay rate 7. When mea- 
surements are performed at time intervals r, we get the effective decay rate 
7eff (■?")■ The full lines represent the survival probabilities and the dotted lines 
their exponential interpolations, according to (j^) . If r = t% < r* one obtains 
QZE. Vice versa, if r — r 2 > t* , one obtains IZE. If r — r* one recovers the 
natural lifetime, according to (||): in this sense, amusingly, t* can be viewed 
as a transition time from Zeno (who argued that a sped arrow, if observed, 
does not move) to Heraclitus (who replied that everything flows), 
opposed Zeno and Parmenides' (Zeno's master) philosophical view 

Sometimes (interestingly) r* does not exist: Eq. (^J) may have no finite 
solutions. In such a case only QZE is possible and no IZE is attainable. This 
is in contrast with some recent claims,CJ according to which the inverse Zeno 
regime is easier to attain than the quantum Zeno one. We can get a qualitative 
idea of the meaning of these statements by looking at the survival probability 
of an unstable system for sufficiently long timeso 

P{t) = \A{t)\ 2 ^Ze-~<\ (11) 



aclitus 
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Figure 3. Z > 1. The full line is the survival probability, the dashed line the renormalized 
exponential e~ 7 ' and the dotted line the asymptotic exponential Ze (a) If P(t) and 
e do not intersect, a finite solution t* does not exist, (b) If P(t) and e intersect, a 
finite solution t* exists. (In this case there are always at least two intersections.) 



where Z, the intersection of the asymptotic exponential with the t = axis, is 
the square modvthriLoJLthe residue of the pole of the propagator (wave function 
renormalizationpo'Ej and will be defined in the next section [Eq. (pi])]. 

A sufficient condition for the existence of a solution r* of Eq. (Q) is that 
Z < 1. This is easily proved by graphical inspection. The case Z < 1 is 
shown in Figure ||(a): P(t) with the property ( |IT|) and e~ 7 * must intersect. 
The other case, Z > 1, is shown in Figure ||: a solution may or may not exist, 
depending on the features of the model investigated. There are also situations 
(e.g., oscillatory systems, whose Poincare time is finite) where 7 and Z cannot 
be defined.c3 The transition from Zeno to inverse Zeno is therefore a complex, 
model-dependent problem, that requires careful investigation. We shall come 
back to this issue in the following sections, where ( [Tl] ) will be derived for a 
particular field theoretical model. 



4 The Lee Hamiltonian 

Some of the most interesting Zeno phenomena, including the transition 
from a Zenja-lo a Heraclitus regime, arise in a quantum field theoretical 
framework.Q'Lnij We will now study the time evolution of a quantum system 
in greater detail, by making use of a quantum field theoretical techniques, and 
discuss the primary role played by the form factors of the interaction. 

Consider a generic Hamiltonian H and an initial normalized state |a). The 
total Hilbert space can always be decomposed into a direct sum TL = Tl a ( S'Hd, 
with Tt a = P a TL and Hd = PdHi where P a — \a)(a\ and Pd = 1 — P a . Let us 
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accordingly split the total Hamiltonian into a free and an interaction part 

H = H a + H int , (12) 

where 

H = P a HP a + P d HP d , H int = P a HP d + P d HP a . (13) 

This decomposition can always be performed,0 even in relativistic quantum 
field theory,E3 the only "problem" being that the decomposition itself depends 
on the initial state \a). Let {\n}} be the eigenbasis of Hq in Ttd 

(a\a) = l, (a\n) = 0, (n\n') = 5 nn ,, (14) 
H \a)=w a \a), H \n) = u> n \n). (15) 

The interaction Hamiltonian H lnt is completely off-diagonal and has nonvan- 
ishing matrix elements only between H a and Tid, namely 

(a\H int \n) = (n\H int \a}* = g n , (a\H int \a) = (n\H int \n') = 0, Vn,n'. 

(16) 



Equations (Il4)-(16) completely determine the free and interaction Hamilto- 



nians in terms of the chosen basis. Indeed we get 

H Q =u a \a)(a\+^2uj n \n)(n\, H int = ^ (g n \a)(n\ + g*\n)(a\) . (17) 

n n 

This is the Lee Hamiltonianil and was originally introduced as a solvable 
quantum field model for the study of the renormalization problem. The in- 
teraction of the normalizable state \a) with the states \n) (the formal sum in 
the above equation usually represents an integral over a continuum of states) 
is responsible for its decay and depends on the form factor g n . 

The Fourier-Laplace transform of the survival amplitude A(t) in (|l|) is 
the expectation value of the resolvent 

G a (E) = -i I dt e' Et A(t) = (al^jjla), 



i 



A(t) = — j dE e^ Et G a (E), (IN) 



the Bromwich path B being a horizontal line IthE =constant> in the half 
plane of convergence of the Fourier-Laplace transform (upper half plane) . By 
performing Dyson's resummation, the propagator reads 

G -< E > - (19 » 
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where the self-energy function 

ME)= ^m^Ml = [ d (20) 

E -lu„ J E-w 

consists only of a second order contribution and is related to the form factor 
g n by the equation 

K °( E ) =^M#i„tN| 2 £(£-^) ^Y,\9n\ 2 5{E-u n ). (21) 

n n 

A comment is now in order. If one is only interested in the survival amplitude 
[or, equivalently, in the expression of the propagator (|l9|)1 and not in the 
details of the interactions g n between \a) and different states \n) with the 
same energy uj n — u>, one can simply replace this set of states with a single, 
representative state \oj), by replacing the Hamiltonian (|l7| ) with the following 
equivalent one 

H = HQ + H int = uj a \a)(a\+ / duj lj\lo){uj\+ \ dw g(u)(\a){w\ + \u){a\), (22) 



where the form factor g(u>) — y K a (u}) and 

|a)(o| + J duj = 1. (23) 

In terms of the Hamiltonian ( p2| ) the self-energy function simply reads 

S4E)= U IfeEafcil! -Up*. (24) 

J hi — uj J hj — lu 



5 Unstable systems 



We consider now the case of an unstable system. The initial state has energy 
uj a > uj g (u) g being the lower bound of the continuous spectrum of the Hamil- 
tonian H) and is therefore embedded in the continuous spectrum of H. If 
— T, a (ujg) < Lu a (which happens for sufficiently smooth form factors and small 
coupling), the resolvent is analytic in the whole complex plane cut along the 
real axis (continuous spectrum of H) . On the other hand, there exists a pole 
.Epoie located just below the branch cut in the second Riemann sheet, solution 
of the equation 

Spolc — U a — Eall(-Epolc) = 0, (25) 
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S a ii being the determination of the self-energy function in the second sheet. 
The pole has a real and imaginary part 

Spoic = uj a + 6u> a - ij/2 (26) 

given by 

Siv a = ReSaiiOEpoic) ~ Re£ (u; + i0+) = p[ dw-^-^-, (27) 

7 = -2ImI] a „(£p io) ~ -2ImE a (w a + iQ+) = 2ng 2 {u a ) 1 (28) 

up to fourth order in the coupling constant. One recognizes the second-order 
energy shift 5uj a and the celebrated Fermi "golden" rule 7JI3 The survival 
amplitude has the general form 

A(t) = A po i e (t) + A cut (t), (29) 

where 

A pole (t) = e -*(«.+*-.)*-r*/a/[ 1 _ V' aU (E pole )}, (30) 

Acut being the branch-cut contribution. At intermediate times, the pole con- 
tribution dominates the evolution and 

P{t)~\A po x e {t)\ 2 =Ze-i\ Z=\l~Z' all (E polc )\-\ (31) 

where Z, the intersection of the asymptotic exponential with the t = axis, 
is the wave function renormalization. We have found ( pd| ) and explicitly de- 
termined Z. 

Notice that, in order to obtain a purely exponential decay, one neglects all 
branch cut and/or other contributions from distant poles and considers only 
the contribution of the dominant pole. In other words, one does not look at 
the rich analytical structure of the propagator and retains only its dominant 
polar singularity. In this case the self-energy function becomes a constant 
(equal to its value at the pole), namely 

G a {E) — ► G™ W (E) = \ = , (32) 

where in the last equality we used the pole equation (^5|). This is the cel- 
ebrated Weisskopf-Wigner approximationtj and yields a purely exponential 
behavior, A(t) = exp(— iE po \ c t), without short- and long-time corrections. 

As is well known, the exponential law is corrected by the cut contribution, 
which is responsible for a quadratic behavior at short times and a power law 
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at long times. In particular, at short times, by plugging (|T^) into ( |l8| ) and 
changing the integration variable rj — Et, Eq. ( |l8| ) becomes 

■Aft) = 7^ / <fy 6 "L ? / v ■ (33) 

The self-energy function ( p4j ) has the following behavior at large energies 

Sa ~ I / = 4<«l^i 2 n tl«> = ^ ^ 00 ( 34 ) 

where we used Eq. (Q) (and assumed the existence of the second moment of 
the Hamiltonian H lnt ). Therefore, the survival amplitude at small times has 
the asymptotic expansion 

A{t) ~^j B dn rf- Ua tn - tVri = 2^J B dr, (r)- t^r, - t m ) ' (35) 



^ = ?±1? +3- (36) 




where 



By closing the Bromwich path in Eq. fl3q ) with a large semicircle in the lower 
half plane, the integral reduces to the sum of the residues at the real poles 
trji t 2 an d the survival probability at small times reads 

P(t) = \A(t)f ~ Hl + Vl-Iwos^^^)] ^ x + = 1 * 

in agreement with the expansion (Q) . Notice that at short times the behavior 
is governed by two "effective" poles which replace the global contribution of 
the cut and the pole on the second sheet. We will come back to this important 
point in the following sections. 

6 Two-pole model and two-pole reduction 

We consider now a particular solvable model: let the form factor in (p2[) be 
Lorentzian 



(38) 

This describes, for instance, an atom-field coupling in a cavity with high finesse 



mirrors E3 (Notice that the Hamiltonian in this case is not lower bounded and 
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-2A to a 2 A 



Figure 4. (a) Form factor g 2 (uS) and initial state energy u) a . (b) Poles of the propagator in 
the complex _E-plane. 



we expect no deviations from exponential behavior at very large times. E9) In 
this case one easily obtains (for ImE > 0) 

= ETiA> ^ 

whence the propagator 

G ^ E) = iE-Jv + iA)-** (40) 
has two poles in the lower half energy plane (see Fig. ||). Their values are 

Ei = u) a + 8u>a-i~, E 2 = -8u a -i(A-~) , (41) 

where 



SuJ a — - ^f- + ^r- 1 ; — t 



2 1 2 V 2H 



with v 2 = w\ + 4A 2 - A 2 . (42) 



7 = A- r 



(Notice that v 2 can be negative.) The propagator and the survival amplitude 
read 

Ex+iA 1 E 2 +iA 1 



G a (E) 



E\ — E2 E — E\ Ei — E2 E — E2 

1 ~ — — 

E-(uj a + 5iu a ) + + E + 6u a + i(A - 7/2) 
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and 

A(t) = (1 - Tl) e -i("«+S"«)t e -yt/2 + ^ e i5a, i e -(A- 7 /2)i 5 ^ 

respectively, where 

1-11 = Res[G a (Ei)] = --L- = Ua \ 6 ;: + l X l12 ] (45) 
l-?t' a {Ex) uj a + 25uj a + i(A - 7) 

is the residue of the pole i?i of the propagator. The survival probability reads 

P(t) = \A(t)\ 2 = Zexp(- 7 i) + 2Re[ft*(l - ft) e -<(<".+2*u.)t] ex p(-Ai) 

+ |^| 2 cxp[-(2A- 7 )t], (46) 

where Z = |1 — 7?.| 2 is the wave function renormalization 

= K + fe a ) 2 + (A- 7 /2) 2 

(w + 2<5w Q ) 2 + (A - 7) 2 • 1 U 

All the above formulas are exact. We now analyze some interesting limits of 
the model investigated. 

6.1 Small coupling 

In the weak coupling limit A <C u) a , A, one obtains from Eq. (E2|) 



LO a + A J CJ a — UJ 

7 = 2A^— + 0(A 4 ) = 2VK) + 0(A 4 ). (48) 

Notice that the latter formula is the Fermi Golden Rule and E\ in ( [4l| ) is 
the "dominant" pole. Indeed, the second exponential in Eq. is damped 
very quickly, on a time scale A" 1 much faster than 7 -1 , whence, after a short 
initial quadratic (Zeno) region of duration A -1 , the decay becomes purely 
exponential with decay rate 7. Notice that the corrections are of order A 2 

K= 2 A2 A2 CJa '^ +°( A4 ) ( 49 ) 
w 2 + A 2 w a + iA 

and the Zeno time is Tz = A -1 ~3> A -1 , i.e. the initial quadratic (Zeno) region 
is much shorter than the Zeno time: in general, the ZjcpDJtime does not yield 
a correct estimate of the duration of the Zeno region. B€3Ej (Beware of many 
erroneous claims in the literature!) The approximation P(t) ~ 1 — t 2 /r| holds 
for times t < A -1 <C tz- 
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6.2 Large bandwidth 

In the limit of large bandwidth A ^> ui a ,A, from Eq. ( ff2| ) one gets 7 = 
2A 2 /A + 0(A -2 ) and in order to have a non trivial result with a finite decay 
rate, we let 

A 2 7 

A — > 00, A — ► 00, with — = — = const. (50) 

A 2 y ' 



In this limit the continuum has a flat band, g(ui) — ^/j/2tt =const, and we 
expect to recover a purely exponential decay. Indeed, in this case one gets 
1Z = and 5uj a — 0, whence 

G a {E) = ^— — , (51) 

E - uj a + «7/2 

so that the survival amplitude and probability read 

A(t) = exp (-iuj a t - ^tj and P(t) = exp(-jt). (52) 

In this case the propagator (|^) has only a simple pole and the survival prob- 
ability ( |52| ) is purely exponential. 

6.3 Narrow bandwidth 

In the limit of narrow bandwidth A<w a ,A, the form factor becomes 

g 2 (uj) = \H{uj) (53) 

and the continuum is "concentrated" in uo = 0. Therefore the continuum as 
a whole behaves as another discrete level and one obtains Rabi oscillations 
between the initial state \a) in ( |2^ ) and this "collective" level (of energy 
uj = 0). In fact one gets 

7 = 0, **_-!& + 0. (54) 

where 



fi=\/A 2 + ^ (55) 



i 
4 

is the usual Rabi frequency of a two-level system with energy difference u) a {— 
Lo a — 0) and coupling A. By (^) the survival amplitude and probability read 

^) = 5( 1 + i) e - ,( - +n " + K 1 -i) e " i! "-° ) ' 

P<0 _!-**>(£). (56) 
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If uj a — 0, the survival probability (56) oscillates between 1 and 0. On the 
other hand, if to a ^ the initial state never "decays" completely. 

Incidentally, notice that the Zeno time is still rz = A -1 and yields now 
a good estimate of the duration of the Zeno region. This is, so to say, a 
"coincidence" due to the oscillatory features of the system. 



6.4 Strong coupling 



Another interesting case is that of strong coupling, A ~ A. This is a typi- 
cal case in which the strong coupling provokes violent oscillations before the 
system reaches the asymptotic regime. In the limit A 3> A,cj q , we get 

^A-^ + OtA- 1 ), 7 = -^ + 0(A- 1 ), K= 1 --^±^ + 0(\-% 

(57) 

whence the survival amplitude reads 
A 



A(t) ~ exp 



.w 
-1 — t 
2 



-t 



1 UJ, 

2 + - 



•tA 



4A 



-i\t 



iA 



4A 



which yields fast oscillations of frequency A damped at a rate A <C A. 



(58) 



6.5 Two-pole reduction 

We now show that the two-pole model introduced in this section is the first 
improvement, after the Weisskopf-Wigner pole, in the approximation of a 
generic quantum field model. First note that, according to the Weisskopf- 
Wigner approximation ([32]), an exponential decay is obtained by considering 
a constant self-energy function E a = —17/2, i.e. a resolvent with a single 
pole with negative imaginary part {E\ in Figure On the other hand, as 
we noted in Sec. ||, the initial quadratic behavior of the survival amplitude is 
governed by two effective poles of the resolvent, which ultimately derive from 
the behavior ([m|) of the self-energy function at infinity 

£„(£)~J-> fOT E^w. (59) 

If one wants to capture this short-time behavior while keeping the exponential 
law at later times, and is not interested in the long-time power-law deviations, 
one can proceed in the following way. The requirement for having an expo- 
nential decay, with decay rate 7 for t — > 00, translates into the behavior of 
the self-energy function for E — > 0, namely in the requirement of having a 
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Weisskopf-Wigner constant self-energy function with negative imaginary part 

S«(0) = -ib. (b > 0) (60) 

The simplest form of the self-energy function satisfying both requirements 
(H) and © is 



1 1/r, 



= ^WTifb = ETiM (61) 

By letting rz = 1/ A and 1 /br% — A, this becomes exactly the self-energy func- 
tion of the two-pole model (p9). Therefore the two-pole model is the simplest 
approximation which yields the short time quadratic behavior together with 
the long time exponential one. 

We call the technique outlined in this subsection "two-pole reduction." It 
is useful if one wants to get a first idea of the temporal behavior of a quantum 
field, keeping information on the lifetime (Fermi golden rule), but also on the 
short-time Zeno region. 

Note that the process outlined above can be iterated to find better ap- 
proximations of the real self-energy function E a (£J) by adding other poles 
and/or zeros. But notice also that this approach does not yield the inverse 
power-law tail. Indeed the latter is essentially due to the nonanalytic behavior 
of the self-energy function at the branching point, a feature that cannot be 
captured by an olomorphic function. 

7 Modification of the Lee Hamiltonian 

We now introduce an interesting modification of the Lee Hamiltonian (|22 
which enables us to look at the Zeno region from a different perspective. The 
Hamiltonian (|22| ) describes the decay of a discrete state |a) into a continuum 
of states \bS) with a given form factor g(ui). According to Eqs. (^) and (|22|), 
the Zeno time is related to the integral of the squared form factor by the 
simple relation 

\= f dug 2 (uj). (62) 

On the other hand, for a two-level system {\a), \b)} with Hamiltonian 

H = X(\a)(b\ + \b)(a\), r z = l/A (63) 

the Zeno time rz is just the inverse off-diagonal element A of the Hamiltonian 
[and, of course, this is in agreement with Eq. (|62]), as shown by Eq. (|5^)]. 
One has therefore a simple system in which the Zeno time is manifest in 
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the Hamiltonian itself. We seek now an equivalent decay model, that shares 
with the two-level model this nice property. To this end, let us add a new 
"intermediate" discrete state \b) to the Lee model. Consider then the Rabi 
oscillation A of the two- level system \a), \b) and let the initial state \a) decay 
only through state \b), i.e. couple \b) to a continuum with form factor g b (u>). 
In other words, the Hamiltonian (E2) is substituted by the following one 



H = Lu a \a) (a\ + u>b\b) (b\ + j diu lu\lu) (lo\ 

+A(|o)(6| + |6)(o|)+ / dug b (w)(\b)(u\ + \u>)(b\). (64) 



We require that this Hamiltonian is equivalent to the original one in describ- 
ing the decay of the initial state \a). To this end, notice that the part of 
Hamiltonian describing the decay of state \b) (and neglecting the coupling 
with \a}) is just a Lee Hamiltonian and gives 

Gb ^ = T 1 -^TJFV = l duJ ^ 

E - uj b - T, b (E) J E-lu 

On the other hand, state \a) couples only to state |6) with a coupling A. 
Therefore the evolution of state \a) is just a Rabi oscillation between state \b) 
dressed by the continuum \ui) and state \a), namely 

G a = G° a + G° a \Gh\G a , (66) 

whence 

Therefore, in the modified model, the self-energy function of the initial state 
| a) is nothing but the coupling A 2 times the dressed propagator G b {E) 

E a (E) = X 2 G b (E) = A ' (68) 

E - uj b - l^ b {E) 

Equation ( pq) is the equivalence relation sought. One has to choose the aux- 
iliary form factor g b (u>) in Eq. ( |6^ ) as a function of the original one g(u>), in 
order to satisfy this relation and get an equivalent description of the decay. 
Our interest in this equivalence is due to the asymptotic behavior of formula 

E a (E) ~ = -4?. for E^^ : (69) 

E T-rE 
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(a) (b) 

Figure 5. The decay of state \a) into a Lorentzian continuum (a) is exactly equivalent 
to a Rabi coupling of \a) with a state \b) that in turn exponentially decays into a flat 
continuum \u>) (b). 



which displays the relation between the coupling A and the Zeno time Tz- 
Thus the Hamiltonian ([54|) explicitly reads 

H = u a \a){a\+u b \b){b\ + — (\a)(b\ + \b)(a\) 

+ J dw uj\uj){lu\ + J ' dw g b {u){\b){w\ + \u)(b\). (70) 

In the equivalent model, therefore, the initial quadratic behavior is singled out 
from the remaining part of the decay: the Zeno region, i.e. the first oscillation, 
is nothing but the initial unperturbed Rabi oscillation between states \a) and 
| b) (which initially "represents" the original continuum as a whole). After the 
initial stage of the decay, the coupling g b {uj) between \b) and \uj) (namely the 
details of the original continuum) comes into play and modifies the initial Rabi 
oscillation with its characteristic time scale. This explains from a different 
perspective the difference, already stressed in previous sections, between the 
Zeno time and the duration of the initial quadratic region. 

As an example, we recover the self-energy function ([39|) of the two-pole 
model, by requiring that be constant 

EftCE) = -uj b -iA , (71) 

which implies 

g b (uj) = \ZA/n and uj b = 0. (72) 

In other words, the auxiliary state |6) is placed at the mean energy of the 
original continuum g(u)) and decays into a flat-band continuum with decay 
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rate 7& = 2A: the decay into a Lorentzian continuum [Fig. |f](a)] is exactly 
equivalent to a Rabi coupling with a level that in turn exponentially decays 
into a flat continuum [Fig. 0(b)]: 



H = uj a \a)(a\ + — (\a)(b\ + \b)(a\) + I du lj\w){u)\ 

±J du(\b)(u\ + \u){b\). (73) 

Moreover, if one restricts one's attention to the subspace spanned by {\a), \ b)}, 
it is easy to showc3 that this Hermitian Hamiltonian reduces to the effective 
non-Hcrmitian one 



H = u a \a){a\ + A (\a)(b\ + |6)(a|) - »A|6)(6| 

= - A ) = - (i + ^I+IW) ) (?4) 

Therefore, if one is interested only in the decay of the initial state \a), the 
study of the two-pole model reduces to the study of this simple non- Hermitian 
2-dimensional matrix. One can reexamine all the results of previous sections, 
just by looking at this matrix. We will not elaborate on this here. 

A final comment is now in order. One can draw a clear picture of the 
two-pole reduction, discussed in the previous section, just by looking at the 
construction of the equivalent model. The first approximation of a real decay, 
the Weisskopf-Wigner approximation, is represented by the simple exponen- 
tial decay of level \b) with its time scale 7 b " 1 . The two-pole approximation 
superimposes an oscillating dynamics with time scale A -1 to the latter, yield- 
ing the initial Zeno region. By complicating the model with the addition of 
other dynamical elements with their characteristic scales, one can construct a 
better approximation of the real decay law. 



8 Zeno Heraclitus transition 



We will now study the Zeno-inverse Zeno transition in greater detail, by 
making use of a quantum field theoretical framework, and discuss the primary 
role played by the form factors of the interaction. The reader should refer to 
the discussion of Sees. |||^ where we introduced the effective decay rate 



7eff(r) 



i 



log P(r) 



■log|.A(T) 



-Re 



log^(r) 



(75) 
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which is a linear function of r for sufficiently small values of r (inside the 
Zeno region) 

7cff(r)~^, for r<l/A, (76) 
T z 

and becomes, with excellent approximation, a constant equal to the natural 
decay rate at intermediate times 

7 off (r) ~ 7 for r > 1/A. (77) 

The transition between Zeno and Heraclitus occurs at the geometrical inter- 
section t* between the curves P(t) and e~ 7 *, solution of the equation 

7eff(r*)=7, (78) 

as shown in Fig. ||. 

Let us corroborate these general findings by considering for example the 
two-pole model studied in detail in Sec. ^|, whose survival amplitude is given 
by Eq. © 

A(t \ = + Sujg + i(A - 7/g) e -i( ttf . +Jtt>tt )t e - 7 t/a 
uj a + 25uo a + i(A - 7) 

+ 5 " a ~ 11 '^ e ^ a t e ^A-,/2)t (?9) 

UJ a + 2Suj a + l(A - 7) 

with 5u> a and 7 given by Eq. (^|). By plugging ([79]) into ^7^ ) one obtains the 
effective decay rate, whose behavior is displayed in Fig. |6j for different values 
of the ratio |w |/A. These curves show that for large values of \ui a \ (in units 
A) there is indeed a transition from a Zeno to an inverse Zeno (Heraclitus) 
behavior: such a transition occurs at r = r*, solution of Eq. (|7§|). However, 
for small values of |w a |, such a solution ceases to exist. 

The determination of the critical value of |oj a | for which the Zeno-inverse 
Zeno transition ceases to take place discloses an interesting aspect of this 
issue. The problem can be discussed in general, but for the sake of simplicity 
we consider the weak coupling limit (small A) considered in Eqs. p^)-(]49|). 
According to the geometrical theorem proved in Sec. |3|, a sufficient condition 
for the system to exhibit an Zeno-inverse Zeno transition is that the wave 
function renormalization Z < 1 . In our case, by making use of Eq. (^9|), this 
condition reads 



\l-K\ 2 = l-2^— Re 
uii + A^ 



u>o — iA 
lu„ + iA 



+ 0{V 
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7cff(»/7 




Figure 6. Effective decay rate 7 e ffM for the two-pole model (^9|), for A = 0.1 and different 
values of the ratio |w a |/A (indicated). The horizontal line shows the "natural" decay rate 
7: its intersection with 7 e ff (t) yields the solution r* of Eq. (|78|). The asymptotic value of 
all curves is 7, as expected. A Zeno (inverse Zeno) effect is obtained for r < t* (t > r*). 
Notice the presence of a linear region for small values of t and observe that t* does not 
belong to s uch linear region as the ratio |u) a |/A decreases. Above a certain threshold, given 
by Eg. (|8l[) in the weak coupling limit of the model (and in general by the condition Z = 1), 
Eq. (p'SI) has no finite solutions: only a Zeno effect is realizable in such a case. 



namely 

c^>A 2 + 0(A 2 ). (81) 

The meaning of this relation is the following: a sufficient condition to obtain a 
Zeno-inverse Zeno transition is that the energy of the decaying state be placed 
asymmetrically with respect to the peak of the form factor (bandwidth) (see 
Fig. |]). If, on the other hand, u> a ~ (center of the bandwidth), no transition 
time t* exists (see Fig. ^) and only a QZE is possible: this is the case analyzed 
inFig.|(a). 

There is more: Equatio n (f79[ ) yields a time scale. Indeed, from the defi- 
nitions of the quantities in flf2|) one gets 7/2 < A — 7/2, so that the second 
exponential in (79) vanishes more quickly than the first one. (The two time 
scales become comparable only in the strong coupling regime: 7 — * A as 
A — > 00.) If the coupling is weak, since 7 = 0(A 2 ), the second term is very 
rapidly damped so that, after a short initial quadratic region of duration 
A -1 , the decay becomes purely exponential with decay rate 7. For t < 1/A 
(which is, by definition, the duration of the quadratic Zeno region), we can 
use the linear approximation ([76[) . When it applies up to the intersection (i.e., 
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\us a \ > A) one gets 

r* ~ 7 r z 2 . (82) 

When uj a gets closer to the peak of the form factor, the linear approximation 
does not hold and the r.h.s. of the above expression yields only a lower bound 



to the transition time r*. In this case the solution t* of Eq. (78) becomes 
larger than the approximation (|82|), eventually going to infinity when condi- 
tion ( |Sl| ) is no longer valid. In such a case, only a QZE is possible and no IZE 
is attainable. This is in contrast with recent claims.c3 



The conclusions obtained for the two-pole model (79) are of general va- 



lidity. In general, in the Lec Hamiltonian (22), for any g(u>), we assume that 
uj a > Lo g (the lower bound of the continuous spectrum), in order to get an 
unstable system. The matrix elements of the interaction Hamiltonian depend 
of course on the physical model considered. However, for physically relevant 
situations, the interaction smoothly vanishes for small values of u — u> g and 
quickly drops to zero for to > A, a frequency cutoff related to the size of the 
decaying system and the characteristics of the environment. This is true both 
for cavities, as well as for typical EM decay processes in vacuum, where tie 
bandwidth A ~ 10 14 — 10 18 s _1 is given by an inverse characteristic lengthEffij 
(say, of the order of Bohr radius) and is much larger than the natural decay 
rate 7 ~ 10 7 - 10 9 s -1 . 

For roughly bell-shaped form factors all the conclusions drawn for the 
Lorentzian model remain valid. The main role is played by the ratio co ag /A, 
where uj ag = u> a — uj g is the available energy. In general, the asymmetry con- 
dition Lj ag < A is satisfied if the energy iv a of the unstable state is sufficiently 



close to the threshold. In fact, from Eq. ( |62| ) one has 

\ = [ do; 5 2 H= ff 2 (w)A, (83) 
r z J 

where u> is defined by this relation and is of order Ci; max , the energy at which 
g(u>) takes the maximum value. For uj a sufficiently close to the threshold lo 9 
one has g(u>) S> g(u> a ), the time scale 7t| is well within the short-time regime, 
namely 



2^.9 2 K) <<; I 
g 2 (Cu)A A 



7*2 - „2,,-\r « T' ( 84 ) 



where the Fermi golden rule 7 = 2itg 2 (uj a ) has been used, and therefore the 
estimate ( |32"| ) is valid. 

On the other hand, for a system such that uj ag ~ A (or, better, uj a ~ center 
of the bandwidth), r* does not necessarily exist and usually only a Zeno 
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effect can occur. In this context, it is useful and interesting to remember 



that, as shown in Sec. 6.3, the Lorentzian form factor (|3S|) yields, in the limit 
g 2 {io) = X 2 S(uj — u a ), the physics of a two level system. This is also true in 
the general case, for a roughly symmetric form factor, when the bandwidth 
A — * 0. In such a case, if ui a = (energy of the initial state at the center of 
the form factor), the survival probability oscillates between 1 and and only 
a QZE is possible. On the other hand, if ui a ^ (initial state energy strongly 
asymmetric with respect to the form factor of "width" A = 0) the initial state 
never decays completely. By measuring the system, the survival probability 
will vanish exponentially, independently of the strength of observation, whence 
only an IZE is possible. 

If one consider the large bandwidth limit of the two-pole model, which is 
equivalent to a Weisskopf-Wigner approximation, the propagator (|5l]) has only 
a simple pole and the survival probability ( |52| ) is purely exponential. There- 
fore measurements cannot modify the free behavior. Indeed, the conditions 
for the occurrence of Zeno effects are always ascribable to the presence of an 
initial non-exponential behavior of the survival probability, which is caused 
by a propagator with a richer structure than a simple pole in the complex 
energy plane. 



9 Continuous measurements 

Most of our examples dealt with instantaneous measurements, according to 
the Copenhagen prescription. Our starting point was indeed Eq. (|6|). How- 
ever, it is always possible to mimic the effect of pulsed measurements in terms 
of the coupling to a suitable system, performing a ccptixuious measurement 
process. This issue has been discussed in other papersj£3o so let us only give 
here an example. Let us add to (|2^) the following Hamiltonian 

H meas (T) = J dtodcu' (\lo)(l;,lu'\ + \lo,lo')(lu\) + J duj' lo'\lu')(lo'\ : (85) 

as soon as a photon is emitted, it is coupled to another boson of frequency w' 
(notice that the coupling has no form factor) . One can show that the dynamics 
of the Hamiltonian (E2) and (pq), in the relevant subspace, is generated by 



H = u a \a)(a\ + J duj (u - iT /2)\uj)(uj\ + J duj g(u)(\a)(uj\ + \uj)(a\), (86) 

and an effective continuous observation on the system is obtained by increasing 
r. Indeed, it is easy to see that the only effect due to T in Eq. (|8^) is the 
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substitution of E (.E) with E a (E + tT/2) in Eq. @, namely, 

G a (£) = -7; „ ,„ : -p /oV (87) 

E — u} a — T, a {E + iT/2) 

For large values of T, i.e., for a very quick response of the apparatus, the 
self-energy function (134) has the asymptotic behavior 



T\ 2 f 2 

-E + i-^-J — if J duj a 2 ^) = -if^' for r-^ 00 - 

[Notice that r — > oo in (|88|) means r ^> A, the frequency cutoff of the form 
factor.] In this case the propagator ( |87j ) reads 

G a (S)~- -A — -, for T^oo (89) 

E - uj a + ij cS (r)/2 

and the survival probability decays with the effective exponential rate (valid 
for T > A) 



7eff(r) ~ -3=. (90) 

Notice the similarity of this result with (Q): T, the strength of the cou- 
pling to the (continuously) measuring system, plays the same role as r _1 , 
the frpgii Piiny of measurements in the pulsed formulation. This is a general 
result E3'H More to this, we have here a scale for the validity of the linear 
approximation ( |90| ) for 7 e g: the linear term in the asymptotic expansion ( |S8| ) 
approximates well the self-energy function only for values of T that are larger 
than the bandwidth A. For smaller values of V one has to take into account 
the nonlinearities arising from the successive terms in the expansion. 

Note that the flat-band case (^l|), yielding a purely exponential de- 
cay, is also unaffected by a continuous measurement. Indeed in that case 
T, a (E) = —ij/2 is a constant independent of E, whence T, a (E+iT/2) = T, a (E) 
is independent of T. The same happens if one considers the Weisskopf-Wigner 
approximation j32| ) : in this case one neglects the whole structure of the propa- 
gator in the complex energy plane and retains only the dominant pole near the 
real axis. This yields, as we have seen, a self-energy function which does not 
depends on energy and a purely exponential decay (without any deviations), 
that cannot be modified by any observations. 

10 Conclusions 

The form factors of the interaction Hamiltonian play a fundamental role when 
the quantum system is "unstable," not only because of the very formulation 
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of the Fermi golden rule, but also because they may govern the transition 
from a Zeno to an inverse Zeno (Heraclitus) regime. The inverse quantum 
Zeno effect has interesting applications and turns out tq_be relevant also in 
the context of quantum chaos and Anderson localization]^ 

Although the usual formulation of QZE in terms of repeated "pulsed" 
measurements a la von Neumann is a very effective one and motivated quite 
a few theoretical proposals and experiments, we cannot help feeling that the 
use of continuous measurements (coupling with an external apparatus that 
gets entangled with the system) is advantageous. 

Both quantum Zeno and inverse quantum Zeno effects have been experi- 
mentally confirmed, ft is probably time to refrain from academic discussions 
and look for possible applications. 
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